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Summary .  Assume  that  each  random  variable  of  a  sequence 

had  a  density  which  is  a  Polya  frequency  function  of  order 
two.  To  this  sequence  we  apply  the  Jifina  sequential  procedure 
to  determine  a  tolerance  interval.  In  this  paper  we  find  some 
sufficient  conditions  on  the  type  of  trend  permissible  for  this 
sequence  which  enable  us  to  show  that  when  the  Ji3fina  procedure 
is  used  the  sampling  will  stop  sooner  and  the  tolerance  interval 
cover  more  of  the  population  (in  a  stochastic  sense)  than  would 
occur  in  the  case  without  trend. 

Similar  considerations  are  shown  to  hold  when  the  sequences 
of  observations  have  densities  which  have  non-decreasing  hazard 


rates 


1. 


The  Jirina  Procedure 


In  the  interest  of  completeness  let  us  define,  in  terms  suf¬ 
ficient  for  our  use,  the  Jirina  procedure.  Let  be  a 

sequence  of  continuous  independent  real  random  variables  (r.v.'s)  not 
necessarily  identically  distributed.  The  triple  (q,  k,  D)  defines 
a  Jirina  procedure  where  ,  k  are  positive  integers  Md  D  is  a 
sequence  of  functions  determined  in  the  following  manner: 

For  given  n  >  t]  let 


(1.1) 


l.n 


<  X.  < 
2,n 


<  X 

n,n 


be  the  ordered  r.v.'s  determined  almost  surely  by  X  , ...,  X  ,  and 

n  In 

for  notational  convenience  write  X  =  (X^,...,  X^) .  Then  for  each 
integer  n  >  q  we  have  D  defined  by 

(1.2)  D(x")  =  U  {x:X.  ^  <  X  <  X.  ] 

1-1, n  -  i,n 


where  the  union  is  over  a  preassigned  set  of  exactly  (n  -  q)  i's  with 
the  proviso  that 

(1.3)  D(X")  C  D(X”'^^)  n  >  q 

We  continue  sampling  until  the  stopping  event 

(1.4)  B(X^'^‘^)  =  [D(X")  =  D(x”'^^)] 
occurs.  In  view  of  this  definition  we  have  also 

Ir  ^ 

B(x"  )  =  n  [X  .  e  D(X")]  . 


(1.5) 


2 


Let  N(X)  be  the  random  sample  size  associated  with  D  which 
we  define  as  the  number  of  observations  n  drawn  when  B(X*^)  occurs 
for  the  first  time.  That  is,  letting  superscript  c  denote  complementation, 

n-1 

(1.6)  CN(x)  =  n]  =  n  b‘^(x^)  n  b(x”)  . 

i=il 

Thus  D(X"<«)  is  the  Jirina  sequential  region  determined  by 
T),  k,  D  .  If  P  is  a  fixed  probability  measure,  the  coverage  of  the' 
region  (with  respect  to  P  )  is  the  r.  v.  on  the  unit  interval  defined 
by 

(1.7)  Q(X)  =  P[D(X^^^^)]  . 

The  case  where  the  identically  distributed  has  been 

studied  previously,  and  it  is  known  that 

k 

(1.8)  PrCQ(X)  >  p]  =  (1  -  p)’l  exp  C-n  Z  pV j  ] 

1 

for  which  approximations  are  known  and  tabulations  have  been  made.  Further 
it  is  known  in  this  case  that 

(1.9)  EN^  =  il  +  kS^  ti>l 

where  S  is  a  constant  which  depends  only  upon  q 

n 

In  this  regard  one  is  referred  to  [3]  and  [4]  for  a  discussion  of  ic-he 
properties  of  (1.8)  and  to  a  translation  of  Jirina's  origineil  paper  in  [2] 
for  its  development . 
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2.  Application  to  Life  Testing 

Suppose  we  are  observing  the  life  lengths  of  successively  produced 
components  with  the  initial  manufacturing  refinements  and  design  im¬ 
provements  being  continually  incorporated  in  their  construction.  We 
assume  the  components  are  being  improved  but  the  degree  of  improvement 
is  not  known  exactly  or  cannot  be  quantified  in  terms  qf  the  life  length 
in  service.  We  also  assume  that  the  degree  of  improvement  of  the  com- • 
ponent  in  the  production  run  will  eventually  reach  a  plateau  of  develop¬ 
ment  beyond  which  it  will  not  progress. 

Suppose  we  apply  the  Jifina  procedure  to  this  sequence  of  observations. 
For  instance,  we  might  agree  to  stop  sampling  when  we  have  obtained  for 
the  first  time  15  observations  which  exceed  the  minimum  life  length  obtained 
in  the  sample.  From  the  tables  and  results  in  [33  we  know  that  if  in  fact 
there  was  no  trend,  i.e.,  the  life  lengths  were  not  being  improved,  then 
we  could  expect  to  stop  in  27.7  observations  and  we  would  have  a  distribution 
of  the  coverage  which  has  the  values  given  in  the  display  below.  On 
the  other  hand,  if  there  is  In  fact  a  very  decided  improving  trend,  it  is 
not  inconceivable  that  we  could  stop  early,  say  less  than  twenty  observations. 
Suppose  for  the  sake  of  comparison  we  consider  the  coverage  of  the 

minimum  of  a  fixed  seimple  of  size  n  =  25  in  the  case  of  no  trend.  Now 
hais  distribution  Pr[Q^  >  P3  =  1  -  p” 


Pr(Q^  >  x) 

Pr(Q  > 

X  =  .8 

.99 

.99 

X  =  .9 

.93 

.97 

X  =  .95 

.72 

.73 

X  =  .99 

.22 

.24 
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It  is  the  purpose  of  this  note  to  set  out  certain  conditions 
under  which  in  the  case  with  trend  the  coverage  is  improved  even  though 
the  sample  size  is  decreased.  The  fact  that  under  certain  conditions 
the  coverage  of  a  sample  with  trend  is  stochastically  larger  than  the 
sample  without  and  the  sample  size  is  stochastically  smaller  for  a 
sample  with  trend  than  a  sample  without  is  the  rationale  for  the  use  of 
the  Jirina  procedure. 

We  now  state  more  precj  ?  ly  our  contention. 

Let  Y^,...,  Y^,...  be  a  sequence  of  independent  non-negative 
r.  v.'s  ,  which  we  may  regard  as  representing  life  lengths,  with  Y^ 
having  continuous  distribution  n  =  1,  2,  . 

We  assume  that 

s  s 

(2.1)  <  ^2  -  •••  ’ 

the  inequalities  meant  in  the  stochastic  sense,  which  is  equivalent  to 

(2,2)  ^2.-  ^2-  •”  • 

© 

Let  X^,...,  X^,...  be  a  sequence  of  independent  r.v.'s  on  the 
positive  real  line  identically  distributed  with  continuous  distribution 
F  ,  and  we  further  aissurae  that 

(2.5)  lim  =  F  . 

n  -•CD 

It  follows  from  (2.2)  and  (2.3)  that  there  exists  a  sequence  of 

s 

order  preseirving  trams  formations  T^,  T^,...  such  that  ^  where 
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we  mean  stochastic  equality  to  a  r.  v.  X  with  distribution  F  ,  that 
is,  ~  ^  further  follows  from  (2.2)  that  for  each  integer  i 

(2.4)  T.  >  T 

1  —  1+1 

We  assume  additionally  that  the  functions  are  differentiable 

and 

(2.5)  T!  <  T'  ^  2 . 

1  —  1+1 

We  make  the  intuitively  appealing  assumption: 

(2.6)  Both  the  X's  and  the  Y's  have  an  increasing  hazard  (failure) 
rate . 

Definition:  A  density  f  with  distribution  F  has  a  (weakly)  increasing 
hazard  rate  (IHR)  iff  6  =  f/l  -F  is  non-decreasing. 

We  do  not  attempt  to  justify  this  assumption  but  it  is  a  common 
one.  We  now  state 

Theorem  1  :  If  the  Jifina  procedure  (l,  k,  D)  is  applied  to  both 
the  X  and  Y  sequences  where  D(X^)  =  [X^  ,00  j  ,  then  the 

random  sample  size  and  the  coverage  with  respect  to  F  , 

say  Q(Y)  ,  for  the  Y  sequence  are,  respectively,  stochastically 
smaller  and  larger  than  the  sample  size  N  and  coverage  with 

A 

respect  to  F  ,  say  Q(X)  ,  associated  with  the  X  sequence. 

The  proof  will  appear  as  a  consequence  of  the  general  results  to  follow. 
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3,  The  Main  Result 

Let  be  independent  random  variables  on  the  real 

line  with  X  having  distribution  F  .  Let  t  be  an  order  pre- 
n  n 

serving  transformation,  i.e.,  continuous  and  strictly  increasing,  and 
for  some  fixed  integer  m  we  define 


Let  us  write  t(X^)  =  (t  X,,.,.,  t  X  )  for  convenience. 

1  n 

Let  D  be  a  Jirina  sequential  tolerance  function  aus  defined  in 
Part  2  with  given  parameters  ii,  k  .  We  make  the  assumption  that 


(3.2  1) 


D(t  x")  =  t [D(X")] 


and 

(3.2.2)  tCD(X^)]  DD(y”)  DD(x") 


We  have 


Lemma  1  :  The  stopping  events  are  invariant  under  t  ,i.e., 
(3.3)  B(tx")  =  B(X'^)  . 

Proof J  This  follows  immediately  from  (2.6)  and  (3.2). 


Lemma  2  ;  Always 

(3.3.1)  B(x”^‘^)  C  B(y”'^^)  n  =  l,  2,...  . 
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Proof:  We  have  three  cases  (i)  n  +  k  <  m  ,  (ii)  n  <  m  <  n  +  k  , 
(iii)  m  <  n  . 


Case  (i) 
Case  (ii) 

Case  (iii) 


By  (3.3)  we  have  B(x"'^*^)  =  B(tx"'^*^)  =  Bd'^'^^)  , 
By  definition  B(x”^*^)  ,=  [D(t  X*')  =  D(t  , 

but  by  (3.2.1)  D(t  X""^^)  D  Dd”"^^)  .  Thus  since 


D(tx'^)  =  D(y”)  ,  we  have  B(X^^^)  which  implies  B(y'^^^)  . 

Now  B(x'''^^)  =  n  CX  .  e  D(x")]  by  (2.7),  but 

j  Ic 

D(y”)  D  D(x’')  by  (3.2.1)  .  Thus  B(x”'*’^)  C  Pi  CX  e  D(y")]  = 

k  .ai  j=l 

=  p  '  ‘>'^>3  • 

3=1 


We  now  have  quite  generally 

Theorem  2  :  Let  N(Y)  and  N(X)  be  the  random  sample  sizes  for  the 
Jirina  procedure  applied  to  the  sequences  Y^,  Y^,...,  Y^  and 
Xu  X2,...  ,  respectively,  defined  above.  Then  N(Y)  is  stochastically 

smaller  than  N(X) 


Proof;  By  (3. 3.1)  we  have  B(x”)  C  B(y”)  .  Hence  taking  complements 

we  thus  have 

n  .  n  . 

(3.3.2)  n  b'^(x^)  d  n  b'^(y^)  , 

i=l  i=l 

where  B(X^)  =  ft)  for  i  <  t)  +  k  .  Thus  we  have 
[N(X)  >  n]  CN(Y)  >  n]  ,  which  proves  the  result. 

We  now  must  discuss  more  specifically  conditions  under  which  we 
can  have  Q(X)  stochastically  smeiller  than  Q(Y) 
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Definition:  A  non-negative  function  f  is  a  Polya  frequency  of  order  2 
(PF^)  iff  f(x)  =  e  where  i|r  is  convex. 


We  now  state  without  proof  some  well  known  facts. 


Fact 

1  : 

f 

is 

IHR 

iff  1  - 

F 

is 

PF^  iff 

F  is 

PF 

^^2  • 

Fact 

2  : 

f 

is 

PFr 

implies 

f 

is 

IHR  . 

Fact 

3  : 

F 

is 

IHR 

iff  Fj^ 

is 

IHR  where 

V"’  = 

1  -  F(-x)  . 

These  statements  are  for  example  proved  by  Barlow,  Marshall  and  Proschan 
in  [1]  . 


We  now  establish 


Lemma  Ji  i  If  H  ,  G  are  non-decreasing  non-negative  functions  and  F 
is  an  IHR  distribution  and  t  is  a  differentiable  function  such  that 
for  every  x 

(5.3.3)  t(x)  >  X  t'(x)  >  1 


Then  for  each  f. 


(3.4) 


t  ^(^t) 


H(y)  G(ty)dF(y)  <  J  H(y)  G(y)  dF(y)  . 


Proof:  Since  G  can  be  approximated  by  an  increasing  sequence  of  linear 

combinations  of  increasing  non-negative  step  functions,  it  is 
sufficient  to  prove  (5.4)  with  G  replaced  by  c(.,a)  where 
c(x,y)=l  if  x>y  and  zero  elsewhere.  Thus  we  need  prove  that 

4'^(?) 

(3.4.1)  J  ,  H(y)  dF(y)  <  J  H(y)  dF(y)  . 

t  (a)  a 
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o 


But  similarly  it  is  sufficient  to  replace  H  by  c(.,b) 

We  then  have  two  cases  (3.4.2)  b  >  a  and  (3.4.3)  t  ^(a)<b<a  , 
So  we  need  prove 


(3.4.2) 


t  I 

J  dF(y)  <  I  ciF(y) 

b  b 


-1. 


which  is  obvious  since  t  (i^)  <  ^  .We  now  prove 


(3.4.3) 


t 


J  dF(y)  <  f  ciF(y) 


which  will  be  done  since  F(b)  <  F(a)  ,  ifawe  show 

-1  -1 

dF(x)  =  Ft  -  Ft  ^(a)  <  J  dF(x)  =  F(5)  -  F(a)  . 

t  "^(a)  a 

Hence  it  is  sufficient  to  prove  that  Ft  ^(x)  -  F(x)  is  a  non¬ 
increasing  function,  but  setting  A(x)  =  -ln[l  -  F(x)J 

Ft  ^(x)  -  F(x)  =  Cl  -  F(x)]Cl  -  exp  {-At  ^(x)  +  A(x)  3  ] 

Thus  we  need  only  establish  that  A(x)  -  A t(x)  is  non-increasing, 
but  6  =  A'  and  6(x)  <  6t(x).t'(x)  ,  which  proves  the  result 

since  6  is  non-decreasing. 


Corollary:  If  H,  G  are  non-increasing  non-negative  functions  and  F  is 

IHR  and  t  is  a  differentiable  function  such  that 

^  t(x)<x  t'(x)>l'  for  each  x  , 

then  for  each  * 


®  CD 

,r  H(y)  G(ty)dF(y)  <  [  H(y)  G(y)  dF(y)  . 

=  I 


(3.5) 
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Proof:  In  (3.4)  replace  y  by  -y  ,  replace  t(x)  by  - 1(- x)  ,  and 
replace  by  -  5  ,  and  the  result  follows  by  a  change  of 
designation  utilizing  fact  3  . 

We  now  prove 

Lemma  4  :  If  , . . , ,  X  are  r.v.’s  with  densities  f ,  ,  . . .  ,  f  which 
are  PF^  then  both  g  and  h  are  PF^  where 

(3.6)  g(x)  =  Pr[x  <  <  ...  <  X^] 

(3.6.1)  h(x)  =  Pr[X,  <  ...  <  X  <  x]  . 

1  n 

o 

Q  Proof:  By  definition 

(T)  CD  CD 

g(x)  =  J  ...J  J  [1  -  F^(x^)]  '^^2  ^3^^3^  '  ^n^\^  “^n 

X  x^ 

Keeping  in  mind  facts  1  and  2,  »  being  PF^  ,  is  IHR  and  hence 

1  -  F^  is  PF^  .  Thus  since  f^  is  PF^  the  product  (1  -  F  )t2  is 
PF^  and  thus  IHR  .  Therefore  the  integral  underlined  is  PF^  .  This 
argument  repeats . 

To  prove  (3.6.1)  simply  replace  X^  by  ”  X^  and  x  by  -  x  and 
keep  fact  3  in  mind. 

We  are  now  in  a  position  to  prove 

Theorem  3  :  If  X^,  X^,...  is  a  sequence  of  r.v.'s  with  densities 
which  are  PF^  and  Y^,  Y^,...  is  a  sequence  determined  by  some 
transformation  t  as  in  (3.I),  then  sufficient  conditions  that  Q(X) 
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be  stochastically  smaller  than  Q(Y)  are 

(a)  t  is  a  differentiable  function  such  that 

t(x)  >  X  t'(x)  >  1  for  all  x 

and  for  some  positive  integer  ti 

D(X")  =  (-00,  X  ,  )  n  >  11 

*  n-Ti+ljn  —  ' 

or 

(b)  t  is  a  differentiable  function  such  that 

t(x)  <  X  t(x)  >  1  for  all  x 

and  for  some  positive  integer  t) 

DCX*^)  =  (X  ,  go  )  n  >  T,  . 

T| ,  n 

Proof:  One  checks  that  for  both  cases  (a)  and  (b)  we  have  (3.2.1)  and 
(3.2.2)  satisfied.  It  is  necessary  to  show  that 

PrCQ(Y)  <  P]  <  Pr[Q(X)  <  p]  for  all  p  e  (0,1) 
but  by  setting 

S^(Y)  =  CQ(y”‘^^)  <  P  ,  N(Y)  =  n  +  k]  n  >  T) 

it  is  sufficient  to  show  that 

Pr  S  (Y)  <  Pr  S  (X)  for  n  >  n  • 

n  —  n  —  ‘ 


m 


Referring  to  a  result  in  [3]  we  have  that 

S^(Y)  =  [Q(y")  <  P  ,  n  ,  Y^  ^  D(Y^~^)  ,  n  e  D(y”)]  . 

In  words  this  formula  merely  says  that  stopping  at  exactly  (n+k) 
observations  and  having  a  coverage  of  less  than  p  can  be  accomplished 
only  by  not  having  stopped  on  or  before  (n  -  1)  observations,  the  nth 
observation  falls  outside  the  tentative  region  and  the  succeeding  k 
observations  fall  inside  the  region  constructed  after  n  observations, 
which  has  a  coverage  of  less  than  p 

Again  we  consider  three  cases  (i)  n+k<m  (ii)  n<m<n+k 
and  (iii)  m  <  n. 


Case  (i)  :  n  +  k  <  m  .  By  (5.I)  we  have  immediately  that 

n-1  .  1  k 

s  (y)  =  CP(tDx^)  <  p  ,  n  B°(x^),  x  i  d(x  ),  n  x  .  e  d(x")  ]  , 

n  —  .  n  .  .  n+.i 

1=T1  j  =  l 


and  by  (3.2.2)  we  have  S  (Y)  C  S  (X) 

n  n 


Case  (ii)  :  n<m<n+k  ,  We  have 


(3.7)  S^(Y) 


[P(tDx")  <  p  , 


n-1  . 

n  B®(x^),  x^  i  D(x”"^), 


m  n+k 

n  X.eD(x''),  n  X.etD(x")]. 

j=n+l  ^  j=m+l  ^ 


Now  we  must  be  more  specific .  In  what  follows  we  shall  consider 
the  case  (a)  ,  D(X")  =  (  -  <D  ,  jj)  »  every  step  may  be  duplicated 

for  case  (h) . 


1 
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Thus  (3.7)  becomes 


S  (Y)  =  {F(tX 
n 


n-1 


.  )<  p  ,  n  B  (X  )  ,  [X  >  X  , 

n-Ti+l.n  —  ’  .  ’  n  n-Ti.n-l  ’ 

'  ’  i=n  '  ’ 


n+k 

n  [X .  <  t  X 

j=m+l  J 


,  ]  ,  n  [X.  <  X  ,  ]} 

n-n+l,n  ^  n-ii+l,n 


Since  for  S  (Y)  to  have  occurred  we  must  have  had  exactly  one  of 
n 

[X.  <  X.  <  ...  <  X.  ]  =  K/ . V  where  (i)  is  a  certain  one  of  the 

1,  1  (1) 

1  d  n  ^ 

permutations  of  the  indices  (l,...,  n)  .  Thus  let  F  (p)  =  ^ 


Pr 


(3.7.1) 


S  (Y)  =  E  J  PrCX,  <...<X.  <y].Pr[y<X.  <. 

(i)  -oo  1  ^n-ti  ^n-11+2 


n+k  n 

TT  F  (ty).  TT  F  (y)  dF 

j=m+l  ^  j=n+l  ^n-Ti+1 


(y) 


but  (3.7.1)  of  the  general  form 

t  ^5) 

J  H(y)  G(ty)  .  f(y)  dy  , 


-  CD 


where  H  and  G  are  monotone  increasing  and  f  is  PF2  so  that  lemma  3 
applies,  and  this  proves  the  result. 

We  remark  that  if  h  =  1  the  second  probability  statement  under  the 
integral  in  (3.7.1)  is  gone,  in  which  case  lemma  3  applies  with  only  the 
assumption  of  IHR  distributions  instead  of  PF^  densities. 

Case  (iii)  :  m  <  n  .  We  have  by  definition 

n-1  .  ,  k 

S  (Y)  =  {p[D(t'^)]  <  P,  n  b'^(y''),  X  eD(Y"~-^),  H  X  e  D(Y*')} 

-i-.,  -i-i  “■'■J 
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but  utilizing  (3.3.2),  and  by  (3.2.2)  that  D(y"~^)  D  D(x'^  ^)  ,  there 
follows 

[X  I  D(y''"^)]  C  [X  e  D(X""^)]  . 

n  n 

Since  we  also  have 

t  D(x")  3  D(Y^)  , 

we  obtain 

n-1  k 

S  (Y)  C  {  PCt  D^X*^)]  <  P  ,  n  b'^(X^)  ,  X  4  D(x”)  ,  n  X  .  e  t  D(X^)  }  , 
n  —  .  n  ’  .  ,  n+n 

i=Ti  j=l 

but  by  case  (ii)  we  see  that  the  probability  of  the  right  hand  side  does 

not  exceed  Pr[S  (X)]  .  This  concludes  the  proof, 

n 

Theorem  4  ;  If  X^,  X^,...  is  a  sequence  of  independent  identically 
distributed  r.v.'s  with  PF^  densities,  and  T^,  T^,...  is  a 
sequence  of  order  preserving  differentiable  transformations  such 
that 


(3.8) 


<  T^’  <  ... 


and  either  (3.8.1)  or  (5.8.2)  holds  where 


(3.8.1) 


T  .  <  T  .  , 

1  —  1+1 


T  .  (x)  -•  X 
1 


(3.8.2) 


T,  >  T.  , 
1  —  1+1 


T.  (x)  -  X 
1 


Then 


?!  =  .  ^2.  " 


(3.8.3) 
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is  a  sequence  of  r.v.’s  for  which  in  case  (3. 8.1)  the  Jirina 
upper  tolerance  interval  (case  (3.8.2)  the  lower  tolerance  interval) 
can  be  established  with  both  stochastically  smaller  sample  size 
and  stochastically  larger  coverage  than  for  the  X-sequence, 

Theorem  5  •  In  theorem  4  the  assumption  of  PF^  densities  can  be  re¬ 
laxed  to  IHR  densities,  if  the  tolerance  intervals  are  restricted 
to  the  minimum  and  maiximum  of  the  sample  in  case  (3. 8.1)  and  (3.8.2), 
respectively . 


Proof:  Let  t^^,  t^,...  be  a  sequence  of  order  preserving  transformations 

such  that 


lim 

m  -•  CD 


t 

m 


t 


m-1 


Since  the  t.'s  are  order  preserving,  they  must  have  non-negative 
J 

derivatives,  and  the  result  is  proved.  But  further  (3.8.1)  is 
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equivalent  with  which  is  equivalent  with 

t  (x)  >  X  by  (3.8.5).  Similarly  (3.8.2)  iff  t.(x)  <  x  . 
^  J 

Thus  we  may  apply  theorem  3  after  applying  each  one  of  the 
t^'s  ,  and  since  at  each  step  the  distribution  of  coverage 
and  sample  size  are  monotone  increasing  and  decreasing, 
respectively,  this  proves  the  result. 


.17 
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